Quantum criticality out of equilibrium in the pseudogap Kondo model 
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We theoretically investigate the non-equilibrium quantum phase transition in a generic setup: the 
pseudogap Kondo model where a quantum dot couples to two-left (L) and right (R)-voltage-biased 
fermionic leads with power-law density of states (DOS) with respect to their Fermi levels jj-l/r, 
Pc,L(R){'^) oc |cj — A'-L{-R)r 1 a-nd < r < 1. In equilibrium (zero bias voltage) and for < r < 1/2, 
with increasing Kondo correlations, in the presence of particle-hole symmetry this model exhibits a 
quantum phase transition from a unscreened local moment (LM) phase to the Kondo phase. Via 
a controlled frequency- dependent renormalization group (RG) approach, we compute analytically 
and numerically the non-equilibrium conductance, conduction electron T-matrix and local spin 
susceptibility at finite bias voltages near criticality. The current-induced decoherence shows distinct 
nonequilibrium scaling, leading to new universal non-equilibrium quantum critical behaviors in the 
above observables. Relevance of our results for the experiments is discussed. 

PACS numbers: 72.15.Qm, 7.23.-b, 03.65.Yz 



Introduction. Quantum phase transitions (QPTs)[lj, 
the continuous phase transitions occur at zero tempera- 
ture due to quantum fluctuations, in strongly correlated 
electron systems have attracted much attention over the 
last three decades. Near the quantum critical points 
(QCPs) associated with QPTs, thermodynamic proper- 
ties exhibit non-Fermi liquid properties and universal 
scalings. Recently, due to high tunability, nano-devices, 
such as: quantum dots in the Kondo regime [3, |3|, of- 
fer a nevif opportunity to study QPTs. In particular, 
understanding QPTs in nano-systems under nonequilib- 
rium conditions has become one of the outstanding emer- 
gent subjects in condensed matter physics with great 
fundamental importance [J-|6|. In Ref. [5|, the authors 
discovered the distinct non-equilibrium profile in trans- 
port near the localized-delocalized QPT of the Kosterlitz- 
Thouless (KT) type in a generic voltage-biased dissi- 
pative resonance-level (quantum dot) from its equilib- 
rium properties at finite temperatures. The current- 
induced decoherence rate smearing out the transition 
shows highly non-linear voltage dependence, resulting in 
these distinct behaviors near QPT. 

In this paper, we investigate the non-equilibrium quan- 
tum criticality in a different class of generic nano-setup- 
the pse udog ap Kondo (PGK) model|7Hll| in a quan- 
tum dot[l4|. We consider a Kondo quantum dot cou- 
ples to two-left (L) and right (R)-fermionic leads with 
a power-law (pseudogap) density-of-states (DOS) which 
vanishes at the Fermi level fJ-L(R) = i^/2, Pc.l(r){'^) °*^ 
oj — ^L{R)Y with < r < 1. Possible realizations of 
the pseudogap leads include: d— wave superconductors 
(r — l)[10j, graphenejlSl (r = 1), one-dimensional Lut- 
tinger systems (r > 0)[8[, and quantum dots embed- 
ded in a Aharonov-Bohm ring (r = 2)[1J|. In equilib- 
rium {V — 0) and for < r < 1/2, with decreasing 
the Kondo couplings the particle-hole (p-h) symmetric 
PGK model exhibits a "true" QPT (distinct from QPT 
of the KT type[15|) from the Kondo screened phase to 



the unscreened local moment (LM) phase[8|, llO|. Near 
QCP separating these two phases, all observables in equi- 
librium exhibit universalpower-law scalings and have 
been extensively studied[10|, lll|. Nevertheless, there is 
still lack of understanding regarding their corresponding 
out-of-equilibrium quantum critical properties. We shall 
address below this issue with a focus on the universal 
nonequilibrium scaling behaviors near QCP. 
The model and the RG approach. The Hamiltonian of 
the particle-hole (p-h) symmetric PGK model reads: 
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where S''"* = fl.r^'afa, S^^, = c^/^fc,^<^/T^'<TCa,fc,<T are 
the spin-1/2 operators of the electron on the dot and in 
the leads, respectively, r are Pauli matrices, and a, a' — 
L/ R, a, a' =^-1 are the lead and spin indices, respectively, 
is the electron creation operator for the lead a with 
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Fermi energies being ^l/r = ±^/2, and fa is the pseud- 
ofermion operator. The conduction electron leads show 
power-law (pseudogap) DOS with respect to their Fermi 
levels /Xi/^, Pc,L(fl)(w) oc Iw - Hl{r)\'' , and < r < 1. 
In the Kondo regime, the single-occupancy constraint of 
the pseudo-fermions is imposed: J2a fa-f^ = 1- Here, 
the dimensionless inter-lead and intra-lead Kondo cou- 
plings are denoted by gLR = NqJlr, and gLL = Qrr = 
NqJll = NqJj^ji, respectively where A^o — jIT ^^'^ ^o 
is the bandwidth cutoff of the leads. For simplicity, we 
consider here the symmetrical Kondo couplings: gap = g. 
In equilibrium, the one-loop RG scaling equation for g 
reads q-j^t^j = rg — 2^^ 10]. The critical Kondo cou- 
pling gc = ^ separates the Kondo (g > gc) from the 
unscreened local moment (LM) phase {g < gc). Much of 
the equilibrium critical properties can be obtained from 
the cutoff dependence of the renormalized Kondo cou- 
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scale being T* = J'ol '^'/"' )^ oc 



with the crossover energy 
5c ^ gl" and the cor- 




FIG. 1: (Color online) Renormalized Kondo coupling g(io) for 
various bare couplings (in units of Do) for r — 0.2 {g^ — 0.1). 
The bias voltage is 1/ = 0.3. Inset: Universal scaling of T/V 
as a function of V/T* with T being the decoherence rate. 



relation length exponent being v = 1/r. At a finite bias 
voltage, however, the chemical potentials (Fermi levels) 
of the two leads are shifted by ±:V/2. Under various RG 
approaches, the Kondo interaction vertices in general de- 
pend not only on the cutoff scale D, but also on the 
electron energy (frequency) [16J, ll8[ . We employ here a 
weak coup ling 1-loop frequency-dependent RG approach 
of Ref. [5|, ll6| which keeps track of energy of the incom- 
ing electrons. For r = our results agree excellently with 
those via a more sophisticated functional RG approach 
in Ref.|19(. Note that our weak coupling theory for the 
p-h symmetric model Eq. [1] works well only for r -^ Q. 
The QCP between LM and Kondo phases disappears for 
r > 1/2, and our theory breaks down for r near 1[10J. 
Note also that the above p-h symmetric QCP is stable 
against p-h asymmetry for < r < r* = 0.375[l3- We 
therefore restrict ourselves to the p-h symmetric model 
for simplicity. The scaling equation for the Kondo cou- 
plings of our model under this approach reads 16|, ll9| : 
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where F is the current-induced decoherence rate ob- 
tained from the imaginary part of the pseudofermion self- 
energy [16|, /" — ^_^^ — is the Fermi function of the a 

lead and ks = h ^ e = 1. Note that in equilibrium at a 
finite temperature T the RG flows of the Kondo couplings 
are cut off by T; while as within the nonequilibrium RG 
approach they are cutoff hy T <^V , a. much lower energy 
scale than T^[16|. Distinct critical behaviors are therefore 
expected [5|. We shah focus below on what these distinct 
nonequilibrium quantum critical behaviors are. 
We first solve Eq.[2]and Eq.[3]for g{uj) self-consistently at 



T = 0. As shown in Fig. 1, for g > {<)gc, the renormal- 
ized Kondo couplings exhibit peaks (dips) at a; = ±V/2, 
indicating Kondo (local moment) phase; while ^(w) is 
completely fiat at criticality g = gc- The qualitative 
nature of these peaks (dips) in (?(w) agree well with 
Ref. [5|, |6J] as signatures of conducting (insulating) behav- 
ior. The height (depth) of the peaks (dips) get shorter 
(shallower) as one reaches to QCP from the Kondo (LM) 
phase. We restrict ourselves to the LM phase {g < gc) 
where the perturbative RG approach is controlled. The 
full analytical solution for g{uj) in the LM phase in the 
limit of 13 —5- is found to be: 



5(w) = g + giiuj) + g2(w), 
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with V = Y^^ uj — Y^^ uq 



Dn 



^T, i" = Y^, and go being 
the bare Kondo coupling. The peaks (dips) of g{uj) near 
cj = ±V/2 shows a power-law behavior: \g{u}) — g{io = 
±V^/2)| oc |w ^ -jl^ with a width of F. We furthermore 
find analytically via Eq. |3]the universal scaling forms for 
g{uj — 0, V), and g{uj = ±V/2, V). These properties will 
be used in the following analysis to determine various 
novel nonequilibrium scaling behaviors in the LM phase: 
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Nonequilibrium decoherence. The current-induced deco- 
herence F which cuts off the RG fiow is the key to under- 
stand nonequilibrium quantum criticality of the model as 
all nonequilibrium observables depend crucially on the 
scaling behavior of F. As shown in Fig.l (Inset), T/V 
in the LM phase exhibits perfect universal V/T* scal- 
ing over a wide range 10~^ < V/T* < 10^. We believe 
this slow crossover which extends over many decades is 
likely related to the large correlation length exponent 
v = 1/r of the model. To gain more insight, we obtain 
analytical approximated form: F/tt « (1 — j)Vg'^{uj = 
t) + j^g'^i^ — 0) where 5(0;) is well approximated by 
a semi-ellipse for —V/2 < ui < V/2 (see excellent agree- 
ment in Fig. 2 (c) between dotted and dashed lines) [5|. 
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where xi^) = J^oo^*^'^*xi.t) with x{t) = 
-(Tc{S'^°*(i) •S'^°*(0)}) are impurity susceptibihties, 

G^ corresponds to the lesser (advanced) component 
of the conduction electron Green's functions with 
constant (DOS) (the effect of the pseudogap leads has 
been taken into account by the renormalized coupling 
giuj)), and T^^,{uj) = T<^,{-uj). The imaginary part of 
the T-matrix at T = is hence given by: 



T^T'^a' (w) = 
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FIG. 2: (Color online) (a). The imaginary part of the T- 
matrix T" (lj) (in unit of ^^) versus V/T* at T = 0. (b). 
h{V) defined in Eq. 9 versus V/T*. Inset: the diagram for 
the T-matrix. (c).The T = nonequilibrium conductance 

o 2 

G{V) (solid lines) normalized to -^^ in the LM phase versus 
V/T* shows distinct scaling from the equilibrium counterpart, 
Q'^i{T _> V') (dashed lines). The dotted line is the analytical 
form via Eq. 11. (d). The scaling of Vxioc versus V/T* with 
Xioc being the local impurity susceptibility. The bare Kondo 
couplings in (a),(b),(c), and (d) are in unit of Do, and r = 0.2. 



Via Eq. [5] the decoherence F at T = is approximated 
as: 



in agreement with the result in Ref. [20| via a different 
perturbative RG approach to the multi-channel Kondo 
model out of equilibrium. For V ~ T = Q, T" {uj)aa' in 
the LM phase exhibits a power-law dip at a; = 0, T {uj) oc 
\uj\^ . For y > 0, this dip is splited into two at w = ±-j 
with the same power-law: T [uj) — T {uj — V/2) oc |w — 
W^. At the dips of r"(w = ±F/2), we find T"(a; = 
V/2) ex (?^(w = 1^/2) shows a distinct nonequilibrium 
scaling behavior as a function of V/T* compared to that 
in equilibrium form T {to = 0) (x g'^ {ui — Q) . To extract 
this different scaling behavior more clearly, we define the 
effective depth of the dips for T"(a; = ±V"/2), estimated 
as: 
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It is clear from Eq. [6] that T /V is an universal scal- 
ing function of V/T* . This well explains the scaling be- 
havior obtained numerically (see Fig. 1 Inset). We ex- 
tract further the asymptotic power-law behaviors of T /V 
as a function of V/T*. For F < F < T*, we have 

^(2^)2'-. For 1/ > F > T*, however, we find 

_ ^ ^[1 „ 2{^)-']. At criticality, F = nglV. The 
scaling behavior of the decoherence F (Fig. 2 and Eq.[S]), 
leading to distinct nonequilibrium scaling behaviors of all 
the observables discussed below, is our central result. 
The conduction electron T-matrix. First, we analyze 
nonequilibrium critical properties from the conduction 
electron T-matrix, defined by G^ ^, ^ — Ga,(T 6a,a' + 

G''S~'T^,.'A^)GTi^)W^ withG«„,,, <(°), be- 
ing the full and bare conduction electron Green's func- 
tion, respectively. The imaginary part of T— matrix 
Im[T{Lo)\ = T {uj) is directly proportional to the experi- 
mentally measurable tunneling density of states (TDOS) 
of our setup. Via renormalized perturbation theory up 
to second order (see Fig. 2 (b) Inset), we have 

T^A^)= E I dng{u)g{uj+n)[x''{mt+X^i^)G^] 



It is clear from Eq. |9]that in the LM phase h{V) follows 
an universal scaling function of V^/T* (see Fig. 2 (b)), and 
has the following asymptotic behaviors: h{V) ~ 1 — (^)'^ 
for F < 1/ <C T*- while for T* < F < F, h{V) w 
2[(^S-)~^ — (w^r)"'^]- The new nonequilibrium scaling 
function h{V) is detectable via STM measurement. 
The nonequilibrium conductance. Next, we turn our at- 
tention to the transport. The nonequilibrium current / 
via the Fermi-Gordon rule readsIS, 161: 



/ = 



T 



I duj\ 
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The current / is computed numerically by Eq. 1101 and 
is approximated asQ: / w ^[{1 - j)Vg^{uj = -j) + 
jVg'^{uj = 0)]. The differential conductance is readily 
obtained numerically via G = ^. In the LM phase, it 
has the analytical approximated form: 



G{V) 



3=L,R 



+ 



37rgg^^ ^^[l + {l + r){^)- 



STT^ff? [1 




16 



< V 
\2T* 



(11) 



(7) As shown in Fig. 2 (c), for T* -^ V ^ Dq, G{V) ap- 



preaches the equihbrium scahng form 
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; while as for V/T* ^ 1 it exhibits a distinct universal 
scaling behavior of V/T* . The perfect scaling behavior 
of G{V/T*) is a direct consequence of the V/T* scal- 
ing in T /V . By contrast, the universal V/T* scaling is 
absent in Ref.[5| as T /V is not a universal function of 
V/T* there. For T < T/ < T*, the conductance behaves 
as: G{V) « H^f^ + b{^f^+'^^"] with a = ^, 

b = ^^f|-(l — f )(%)'^, which shows a pref actor reduc- 
tion in the first term oc V"^^ with respect to its equi- 
librium form G"=«(F ^ T <^ T*) ^ ^[^?'' and a 
sub-leading correction with an anomalous power-law be- 
havior oc l/2''+2'- ^ For y > r » r*, however, we find 
G{V) « ^ [1 - p( 2^)-'- - g(^)-2l with p = 1 - r 



and <Z=f(l — f)(l— ?'){^) 2, which deviates signifi- 
cantly from its equilibrium form G'^'^iV -^ T ^ T*) « 
■^^Yg-[1 — 2 (y^ )"'■]. It is worthwhile emphasizing that due 
to the very different role played by the bias V and tem- 
perature T, G{V) follows a completely different scaling 
function from its equilibrium form G^'^{V -» T) over the 
full range of V/T* (see Eq. [11] and Eq. [TH) though it 
tends to converge with its equilibrium form for V <^T*. 
Local spin susceptibility. We furthermore analyze 
the scaling behaviors of the local spin susceptibility 
XiociV) = ^^\h^o in the LM phase with h being a small 
magnetic field, M being the magnetization M = ^-r^- 
Following Ref. 0, i, [3, for T/ > /i ^ 0, we find M w 

_ The approximated form for Vxioc reads 

(see Fig. 2 (d))^: Vxioc « .,^^=oAT-Si^=v/2y^ ' 
For r <^ 1/ <C T*, xioc exhibits an anomalous power-law 



/ 



du}g (a 



behavior: 



Xlo 



oc 



with 7]-^ 



-- 2r2 
Txioc 



distinct from 
(T < T*) oc 



its equilibrium constant behavior 

(ffc - 5)*^ [ill- For y > r > T*, however, we find the 
local susceptibility acquires a power-law correction to the 
Curie behavior: Vxioc ~ \-V l^xioc and l^xioc oc pr^s^ 
with an anomalous exponent Ar;^ — — r; while as its cor- 
responding equilibrium form shows a different anoma- 
lous power-law behavior: xiodT » T*) oc rpi-„^ with 
Vx ~ ''^ idjlUl- At criticality {g — gc), xioc shows per- 
fect Curie law behavior: xioc ex l/V. These distinct 
nonequilibrium signatures near QCP are detectable in 
local susceptibility measurement. 

Conclusions. In summary, via a controlled frequency- 
dependent renormalization group approach we have in- 
vestigated the quantum phase transition out of equilib- 
rium in the pseudogap Kondo quantum dot. At zero 
temperature and finite bias voltage, we discovered in the 
local moment phase the new quantum critical behaviors 
in the T-matrix, conductance, and local spin suscepti- 



bility that are distinct from those in equilibrium and at 
finite temperatures. The key to explain these differences 
lies in the fact that the current-induced decoherence at a 
finite bias voltage (out of equilibrium) acts quite differ- 
ently from that at a finite temperature but zero bias (in 
equilibrium) , resulting in distinct nonequilibrium behav- 
ior near the quantum phase transition. Our predictions 
open up a new perspective both theoretically and experi- 
mentally in the study of the Kondo dot coupled to exotic 
leads with pseudogap density of states. 
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